WAVELET FILTERING 
Technical Field of the Invention 

The present invention relates generally to wavelet filters and, in particular, to 
wavelet filtering of digital images. 

Background 

The field of digital data compression and in particular digital image compression 
has attracted great interest for some time. 

In the field of digital image compression, many different techniques have been 
utilised. In particular, one popular technique is the JPEG standard, which utilises the 
discrete cosine transform to transform standard size blocks of an image into 
corresponding cosine components. In this respect, the higher frequency cosine 
components are heavily quantised so as to assist in obtaining substantial compression. 
The heavy quantisation is an example of a lossy technique of image compression. The 
JPEG standard also provides for the subsequent lossless compression of the transform 
coefficients. 

Recently, the field of wavelet transforms has gained great attention as an alternative 
form of data compression. The wavelet transform has been found to be highly suitable in 
representing data having discontinuities such as sharp edges. Such discontinuities are 
often present in image data or the like. 

Typically, data compression using wavelet techniques is a two step process. It 
comprises, firstly, a transform phase, during which the wavelet transform of the data set is 
calculated, and secondly a subsequent coding stage during which the resultant data set 
from the transform operation is separated into segments which are then coded using a 
specific coder. In decompression, the reverse occurs, with coded blocks being first 
decoded, and subsequently the inverse wavelet transform being applied to generate the 
final decompressed output. 

Typically, the wavelet transform operation utilises a number of wavelet filters. 
These filters are typically digital filters, such as finite impulse response (FIR filters) 
which ideally are linear and require a small number of taps. These filters process the 
digital images in bands and typically require a large external store for storing overlapping 
pixels, coefficients, and/or intermediate results for processing the next band. The 
bandwidth to the external store is typically limited, which restricts the rate of reading and 
writing of overlapping pixels, coefficients, and/or intermediate results to the store. 
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Disclosure of the Invention 

It is an object of the present invention to substantially overcome, or at least 
ameliorate, one or more disadvantages of existing wavelet filters. 

According to a first aspect of the invention, there is provided a method of wavelet 
5 filtering a digital image, the digital image comprising a plurality of pixels arranged in one 
or more bands of pixels, the method performing the following steps for each one of said 
one or more bands of said pixels: retrieving a plurality of previous partial results from a 
first or second storage; inputting a current group of adjacent said pixels; computing 
coefficients utilising said current group of adjacent pixels and said plurality of previous 

10 partial results; outputting said coefficients; and storing a plurality of current partial results 
in said first storage if said current group is not the last group in the band, in said second 
storage if said group is the last group in the band. 

According to a second aspect of the invention, there is provided an apparatus fox- 
wavelet filtering a digital image, the digital image comprising a plurality of pixels 

15 arranged in one or more bands of pixels, the apparatus comprising: means for retrieving, 
for each one of said one or more bands of pixels, a plurality of previous partial results 
from a first or second storage; means for inputting, for each one of said one or more 
bands of pixels, a current group of adjacent said pixels; means for computing coefficients, 
for each one of said one or more bands of pixels, utilising said current group of adjacent 

20 pixels and said plurality of previous partial results; means for outputting, for each one of 
said one or more bands of pixels, said coefficients; and means for storing a plurality of 
current partial results in said first storage if said current group is not the last group in the 
band, in said second storage if said group is the last group in the band. 

According to a third aspect of the invention, there is provided a computer readable 

25 medium comprising a computer program for wavelet filtering a digital image, the digital 
image comprising a plurality of pixels arranged in one or more bands of pixels, the 
computer program comprising: code for retrieving, for each one of said one or more 
bands of pixels, a plurality of previous partial results from a first or second storage; code 
for inputting, for each one of said one or more bands of pixels, a current group of adjacent 

30 said pixels; code for computing coefficients, for each one of said one or more bands of 
pixels, utilising said current group of adjacent pixels and said plurality of previous partial 
results; code for outputting, for each one of said one or more bands of pixels, said 
coefficients; and code for storing a plurality of current partial results in said first storage 
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if said current group is not the last group in the band, in said second storage if said group 
is the last group in the band. 

According to a fourth aspect of the invention, there is provided a wavelet filter for 
wavelet filtering a digital image, the digital image comprising a plurality of pixels 
5 arranged in one or more bands of pixels, the filter comprising: a first storage for storing a 
plurality of previous partial results; a second storage for storing a plurality of previous 
partial results; a controller for selecting the plurality of previous partial results from said 
first or second storage; a pixel input mechanism for inputting a current group of adjacent 
said pixels; a lifting lattice of multiplier and adder units for computing coefficients 
10 utilising said current group of adjacent pixels and said selected plurality of previous 
partial results; output means for outputting said coefficients; and a controller for storing a 
plurality of current partial results in said first storage if said current group is not the last 
group in the band, in said second storage if said group is the last group in the band. 
Brief Description of the Drawings 
15 A number of preferred embodiments of the present invention will now be described 

with reference to the drawings, in which: 

Fig. 1 is a schematic representation of a K-line 5-3 wavelet filter (with K = 8) for 
filtering a digital image in accordance with a first embodiment; 

Fig. 2 is a schematic representation of the cascading of a wavelet filter of Fig. 1 in 
20 time; and 

Fig. 3 is a flow diagram of a method for use in wavelet filtering in accordance with 
a second embodiment; and 

Fig. 4 is a schematic representation of a general-purpose computer for use in 
implementing the method of Fig. 3. 
25 Detailed Description including Best Mode 

Where reference is made in any one or more of the accompanying drawings to steps 
and/or features, which have the same reference numerals, those steps and/or features have 
for the purposes of this description the same function(s) or operation(s), unless the 
contrary intention appears. 
30 The principles of the preferred embodiments described herein have general 

applicability to the implementation of K-line wavelet filters. However, for ease of 
explanation, the embodiments has been described with reference to a K-line 5-3 wavelet 
filter for use in a forward discrete wavelet transform of digital images. However, it will 



541400US.doc 



be readily evident that the invention is not limited thereto. For instance, the wavelet filter 
can also be used in K-line 9-7 wavelet filters. Furthermore, it would be apparent to a 
person skilled in the art, that the wavelet filter can be modified for use in an inverse 
discrete wavelet transform. Still further, it would be apparent to a person skilled in the art 
that two such wavelet filters can be combined in known manner to provide a 2- 
dimensional discrete wavelet transform or modified for use in a 2-dimensional inverse 
discrete wavelet transform. Moreover, the wavelet filter can be used in many different 
applications. For examples of the many different applications of wavelet analysis to 
signals, reference is made to a survey article entitled "Wavelet Analysis" by Bruce et. al. 
appearing in IEEE Spectrum, October 1996, pages 25 to 26. For a discussion of the 
different applications of wavelets in computer graphics, reference is made to "Wavelets 
for Computer Graphics", page 5, 1. Stollinitz et. al. published 1996 by Morgan Kaufhann 
Publishers, Inc. 

As mentioned above, the embodiments of the invention can be used in a discrete 
wavelet filter transform as well as an inverse discrete wavelet transform. For sake of 
simplicity, the terra wavelet transform used herein or variations thereof, is taken to 
include both a forward and/or an inverse wavelet transform, unless the contrary intention 
appears. Similarly, the term pixel as used herein or variations thereof, is taken to taken to 
refer to the original pixels or wavelet coefficients of a digital image, unless the contrary 
intention appears. Similarly, the term digital image or variations thereof is taken include 
an original image or a sub-band of that original image, or one or more associated sub- 
bands of wavelet transformed coefficients of that original image. 

Turning now to Fig. 1, there is shown a schematic representation of a K-line 5-3 
wavelet filter (with KL = 8) for filtering a digital image in one dimension in accordance 
with a first embodiment. The digital image comprises a plurality of pixels arranged in a 
plurality of columns in the vertical direction and a plurality of rows in the horizontal 
direction. , 

The filter 100 processes one horizontal band of the digital image after another. The 
digital image comprises one or more such horizontal bands, with the horizontal bands 
comprising a plurality of said rows of said pixels. The filter 100 initially commences 
processing the uppermost band of the image, and then processes the next adjacent band of 
the image, and so on until the lowermost band of the image is processed. Within each 
band, the filter 200 processes one column of pixels after another. The filter 100 first 
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processes the left column of the band, and then the next adjacent column of the hand, and 
so on until the right column of the band is reached. For ease of understanding, the pixels 
in each column of a band from the topmost pixel to the lowermost pixel are represented as 

Xo.Xi.X2. --XK--I- 

5 Within each column of a band, the filter 100 processes a first group of adjacent 

pixels (commencing with the topmost group) and then processes a second group of 
adjacent pixels adjacent the first group. 

In this way, the filter 1 00 performs a one dimensional discrete wavelet transform in 
the vertical direction. In another embodiment, the filter 100 can also be adapted to 

10 perform a one dimensional discrete wavelet transform in the horizontal direction by 
processing vertical bands one after another in similar fashion. In still another 
embodiment, the filter can also be adapted to provide a 2-dimensional wavelet transform. 
In the latter case, one filter 100 processes one column at a time first, and a second filter 
100 after that collects the output of the first filter 100 and processes one row at a time. 

15 The pixels ( e.g. x 0 ,X} r X2,X3 - X4 r X5_x^X7 f X8.x^xjo,x/ /i x/2 l x/3 > Xj4andxjs) of a column 

of a band are fed to the filter 1 00. These pixels are fed in parallel to the wavelet filter 1 00 
via a number 2K (with K = 8) of input lines 101. The 16 input lines 101 are connected to 
eight 2-input multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 103-6, 103-7, and 103-8. 
The multiplexers can select either the bottom eight lines or top eight lines for input of the 

20 groups of pixels (x 0l x ft x 2> x 3 _ X4. xs. x& and xi) or (xs. x?, x /0 , x, Ji x !2> x ]3% x I4 and xjs) 
respectively. The multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 103-6, 103-7, and 103- 
8 are connected to a lifting lattice of multiplier and adder units 102, 104, 106, 108, 110, 
1 12, 114, and 1 16, which compute the high pass coefficients and low pass coefficients on 
the top 8 lines (ie. coefficients (xg, xg_ xjo, xu, xj2. *J3. xu and x^)) or the bottom S lines ( 

25 ie. coefficients (xo. Xi. x^ x 4 , x 5 , X6. and x 7 )). The lifting lattice of multiplier and adder 
units 102, 104, 106, 1 OS, 110, 112, 114, and 116, perform the following calculations: 

d n =x i« +a { x 2n-i +x in + i) n = 0,\,...,K-\ (1) 

and, 

30 s n = x 2n+l + p(d n + d nH ) n = -1,0,. - 2 (2) 

where d„ are all the high pass coefficients, s„ all the low pass coefficients, a = -0.5, and (3 
= 0.25. 
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For example, the multiplier and adder units 102, 104 106, 108, 110, 112, 114, and 
116 perform the following calculations on the selected input pixels (x 0 . x L x 2 , x 3 . x 4 , x 5 x 6i 
and x?) to compute the high pass coefficients d n and the low pass coefficients s n : 



unit 102 


d a = x 0 + cz(x_ l +x 1 ) 


(3) 


unit 104 


d } -x 2 + a(x l + x 3 ) 


(4) 


unit 1 06 


d 2 = x 4 + a(x 3 + x s ) 


(5) 


unit 108 


d 3 = x 6 + a(x 5 +x 7 ) 


(6) 


unit 110 


-J-, = + J3{d^ +d 0 ) 


O) 


unit 112 


s 0 =x x + /3{d,+d x ) 


(S) 


unit 114 


=x J +j3(d i +d 2 ) 


(9) 


unit 116 


s i = X s + Pid 2 + d 3 ) 


(10) 



The high pass and low pass coefficients d R and s a are then passed to other discrete 
wavelet transform circuits or entropy coder (not shown in the figure). 

The multiplier and adder units 1 10 and 1 02 also accept inputs from either the partial 
result registers 109 or external buffer 107, depending on the control on the 
multiplexer 105. Namely, the multiplier and adder units 110 and 102 accept input from 
the external buffer 1 07 when computing the first 8 lines and accepts input from the partial 
result registers 109 for the last 8 lines. 

A controller 1 1 1 controls the storing of a last input pixel and a last high pass 
coefficient (e.g. x 7 and d 3 ) to the partial result register 109 or external buffer 107 
depending on which lines are being processed. Namely, the last input pixel and the last 
high pass coefficient of the first 8 lines are stored in the partial result register 109. 
However, when the lifting lattice of multiplier and adder units is computing coefficients 
for the last 8 lines, the last input pixel and last high pass coefficient are stored in the 
externa] buffer 107. The last input pixel and last high pass coefficient are required for 
computation of the low and high pass coefficients of the next adjacent group of pixels 
within the column. 

For example, in the case where the filter 100 is currently processing the pixels (x$. 
xp, x,o, x u . x n . xj 3 . xi4 . and x ls ), the multiplier and adder units 1 10 and 102 accept as input 
the last pixel and coefficient x 7 and d 3 from the partial result registers 109, which were 

541400US.doc 



previously stored in the partial result register 109 during the previous computation of the 
pixels (xo. xi_ x 2 , xj, x 4 , x 5: xg. and X7) of the same column and band. In the case where the 
filter 100 is currently processing the pixels (xo to x?), the multiplier and adder units 110 
and 102 accept as input the last pixel and coefficient x/ s and d 7 from the external 
5 buffer 107. These inputs xi 5 and d.7 were previously stored in the external buffer 107 
during the previous computation of the pixels (x s to xjs) of the same column but in the last 
previous band. In the case where there are no previous pixels (e.g. topmost band) the 
boundary is symmetrically reflected. 

Figure 2 illustrates the cascading of the wavelet filter of Fig. 1 in time. In cycle 0, 

10 for the first set of inputs the wavelet filter 100 either symmetrically reflects the pixels 
near the boundary or reads the required partial results from external buffer 107. The 
filter 100 then stores the partial results into the partial results register 109, and outputs the 
low and high pass coefficients. In cycle 1, the lifting lattice is 'moved' to the next set of 
input lines, reads input from partial results register 109, and outputs the next low and high 

15 pass coefficients. The partial results are stored in the partial results register 109 and 
external buffer 107. 

The first embodiment has been described with respect to a filter 100 having two 
groups of eight input lines 101 and eight two-input multiplexers 103-1 to 103-8. It would 
be apparent to a person skilled in the art, that the number of groups can be increased (e.g. 

20 N = 4 groups) with a corresponding modification of the eight multiplexers, (ie. eight N- 
input multiplexers). The external buffer 107 is a slow memory storage and requires at 
least a memory size equivalent to the size of an input pixel and a high pass coefficient 
times the number of columns. The filter 1 00 stores in external memory the last input 
pixel and the last high pass coefficient xj s and d?of each column of each band for further 

25 processing in the next band. On the other hand, the partial results register is a fast local 
storage and requires only a memory size of an input pixel and high pass coefficient. 

The effective value of K' (the width of the band) can be increased by increasing the 
number of groups of pixels (eg. two groups, K' = 2xK — 16) or (four groups K' = 4xK 
=32)). In the present embodiments, this has the advantage of reducing the number of 

30 accesses to the slow external memory. Namely, where there are multiple groups, the last 
pixel and last high pass coefficient can be stored in the fast local memory, namely the 
partial results register 109, rather than the slow external memory 107, thus enabling the 
wavelet transform to be performed more efficiently. Of course, the last pixel and last 
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high pass coefficient of the last group is required to be stored in external memory. 
Furthermore, the amount of hardware required to implement the filter is not significantly 
increased. Still further, the filter has the capability of filtering with different values of K\ 
In an alternate embodiment, the multiplexers 103-1, 103-2, 103-3, 103-4, 103-5, 
5 103-6, 103-7, and 103-8 may be dispensed with. In this case, the multiplexing is done by 
pixel fetching mechanism implicitly. For instance, the pixels are fed to the lifting lattice 
from a frame buffer by an address generator. 

Turning now to Fig. 3, there is shown a method for use in wavelet filtering in one 
dimension in accordance with a second embodiment. The method comprises in part a 

10 sub-procedure 300 for use in the wavelet filtering of a digital image. As mentioned 
previously, the digital image comprises a plurality of pixels arranged in a plurality of 
columns in the vertical direction and a plurality of rows in the horizontal direction. The 
digital image comprises one or more horizontal bands, each of which comprise a plurality 
of said rows of pixels. The method processes one horizontal band of the digital image 

15 after another. The method initially commences processing the uppermost band of the 
image, and then processes the next adjacent band of the image, and so on until the 
lowermost band of the image is processed. Within each band, the method processes one 
column of pixels after another. The method first processes the left column of the band, 
and then the next adjacent column of the band, and so on until the right column of the 

20 band is reached For ease of understanding, the pixels in each column of a band from the 
topmost pixel to the lowermost pixel are represented as xo. x J: x 2 , ... Xk--i- 

The method of wavelet filtering calls the sub-procedure 300 for processing each 
column of a band. Within each column of a band, the sub-procedure 300 processes in 
turn a first group of K adjacent pixels, a second group of K adjacent pixels adjacent the 

25 first group , a third group of K adjacent pixels adjacent the second group, and so on to the 
last group of pixels. The top and bottom pixel of the first and last group of pixels 
respectively, being located adjacent the boundary of the band. 

In this way, the method of wavelet filtering performs a one dimensional discrete 
wavelet transform in the vertical direction. In another embodiment, the method can also 

30 be adapted to perform a one dimensional discrete wavelet transform in the horizontal 
direction by processing vertical bands one after another in similar fashion. In still another 
embodiment, the method can also be adapted to provide a 2-dirnensional wavelet 
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transform. In the latter case, the method processes one column at a time first, and the 
method after that collects the output and processes one row at a time. 

The sub-procedure commences at step 301 and proceeds to step 302 where a 
counter M is set to one. As mentioned previously, the sub-procedure processes the pixels 
within a column of a band. The sub-procedure continues to a loop 306, 308, 310, 312, 
314, and 316, where groups of a number K adjacent pixels are retrieved and processed 
one after another. A first group of K adjacent pixels (nearest the top boundary) is 
processed during the first pass of the loop 304 to 313, a next group of K adjacent pixels 
(adjacent the first group) is processed during the second pass of the loop 304 to 313, and 
so on The processing of one group of pixels after another is achieved by incrementing 
the counter M. The sub-procedure terminates once K' pixels have been processed, when 
M.KL = K\ Preferably, K = 4 and K'= 4 or S. 

After the counter M has been set to one, the sub-procedure proceeds to step 304. hi 
step 304, the method reads the corresponding previous partial results from an external 
buffer. These previous partial results having being stored in the external buffer during the 
processing of the sub-procedure of the same column in the last previous band of pixels. 
In the case where there are no previous pixels (e.g. topmost band) the boundary is 
symmetrically reflected. 

After step 304, the method proceeds to loop 304 to 316. This loop will be 
described with respect to the general case where M = M\ During step 306, the sub- 
procedure reads the M'th group of K adjacent pixels and then computes the low and high 
pass coefficients in step 308 in accordance with the equations (1) and (2) described 
above. In computing these high and low pass coefficients the method takes as input the 
M'th group of K adjacent pixels. It also accepts input from either the external buffer or a 
local register which stores the last pixel and the last high pass coefficient. The step 308 
outputs to other discrete wavelet transform or entropy encoding procedures (not shown in 
the figure) . 

For example, in the case where the sub-procedure is currently processing the pixels 
(X8. *s>. xjo, x n , x I2 , x !3t X14 , and x JS ), the sub-procedure step 308 accepts as input the last 
pixel and high pass coefficient x 7 and d 3 from the partial result registers, which were 
previously stored in the partial result register during the previous computation of the 
pixels (x 0 . xj f x 2 , xj, x 4 , x 5r x 6 . and x 7 ) of the same column and band, hi the case where the 
sub-procedure is currently processing the pixels (x 0 to x 7 ), the sub-procedure step 308 
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accepts as input the last pixel and high pass coefficient from the external buffer. These 
latter inputs were previously stored in the external buffer during the previous computation 
of the last group of pixels of the same column but in the last previous band. In the case 
where there are no previous pixels (e.g. topmost band) the boundary is symmetrically 
5 reflected. 

Ater step 308, the method proceeds to a decision block 310, where a check is made 
whether M.K = K\ If the decision block returns false, then the method proceeds to 
step 312, where the last pixel and last high pass coefficient are stored in the local partial 
results register for use in the computation 308 of the next group of pixels. 

10 In the next step 314, the counter M is incremented and the last pixel and the last 

high pass coefficient are read 316. The method then returns to step 306 for the next pass 
of the loop 306 to 316. In the event, the decision block 310 returns true, then the last 
pixel and the last high pass coefficient are stored in the external buffer and the sub- 
procedure finishes for the current column of the band. 

15 As can be seen, the loop can be repeated as many times as desired to increase the 

value of K\ 

The aforementioned preferred method comprise a particular control flow. There are 
many other variants of the preferred method which use different control flows without 
departing the spirit or scope of the invention. Furthermore one or more of the steps of the 

20 preferred method may be performed in parallel rather sequential. 

The method of Fig. 3 is preferably practiced in dedicated hardware, such as shown 
in Fig. 1. However, the method of Fig. 3 can also be practiced using a conventional 
general-purpose computer system 400, such as that shown in Fig. 4 wherein the processes 
of Fig. 3 may be implemented as software code, such as an application program executing 

25 within the computer system 400. In particular, the steps of method of Fig. 3 are effected 
by code instructions in the software that are carried out by the computer. This method 
has particular application in computer systems restricted to large slow external buffers 
and small fast local memories, 

The software may be stored in a computer readable medium, including the storage 

30 devices described below, for example. The software is loaded into the computer from the 
computer readable medium, and then executed by the computer 
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The use of the computer program readable medium together with the software in 
the computer system preferably effects an advantageous apparatus in accordance with the 
embodiments of the invention. 

The computer system 400 comprises a computer module 401, input devices such as 
5 a keyboard 402 and mouse 403, output devices including a printer 415 and a display 
device 414. A Modulator-Demodulator (Modem) transceiver device 416 is used by the 
computer module 401 for communicating to and from a communications network 420, for 
example connectable via a telephone line 421 or other functional medium. The 
modem 416 can be used to obtain access to the Internet, and other network systems, such 

10 as a Local Area Network (LAN) or a Wide Area Network (WAN). 

The computer module 401 typically includes at least one processor unit 405, a 
memory unit 406, for example formed from semiconductor random access memory 
(RAM) and read only memory (ROM), input/output (I/O) interfaces including a video 
interface 407, and an I/O interface 413 for the keyboard 402 and mouse 403 and 

15 optionally a joystick (not illustrated), and an interface 408 for the modem 416. A storage 
device 409 is provided and typically includes a hard disk drive 410 and a floppy disk 
drive 411. A magnetic tape drive (not illustrated) may also be used. A CD-ROM 
drive 412 is typically provided as a non- volatile source of data. The components 405 
to 413 of the computer module 401, typically communicate via an interconnected bus 404 

20 and in a manner which results in a conventional mode of operation of the computer 
system 400 known to those in the relevant art. Examples of computers on which the 
embodiments can be practised include IBM-PC's and compatibles, Sun Sparcstations or 
alike computer systems evolved therefrom. 

Typically, the application program of the preferred embodiment is resident on the 

25 hard disk drive 410 and read and controlled in its execution by the processor 405. 
Intermediate storage of the program and any data fetched from the network 420 may be 
accomplished using the semiconductor memory 406, possibly in concert with the hard 
disk drive 410. In some instances, the application program may be supplied to the user 
encoded on a CD-ROM or floppy disk and read via the corresponding drive 412 or 41 1, 

30 or alternatively may be read by the user from the network 420 via the modem device 4 1 6. 
Still further, the software can also be loaded into the computer system 400 from other 
computer readable medium including magnetic tape, a ROM or integrated circuit, a 
magneto -optical disk, a radio or infra-red transmission channel between the computer 
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roodule 401 and another device, a computer readable card such as a PCMCIA card, and 
the Internet and Intranets including email transmissions and information recorded on 
websites and the like. The foregoing is merely exemplary of relevant computer readable 
mediums. Other computer readable mediums may be practiced without departing from 
the scope and spirit of the invention. 

The aforementioned embodiment shown in Fig. 1 can be modified for use in a 5/3 
inverse transform wavelet filter. In this modified embodiment, the 5/3 inverse transform 
wavelet filter calculates the original pixels from the low pass coefficients and high pass 
coefficients in accordance with the following formulae: 

x 2 „=s„- p{d n _,+d n ) n = 0,l,...,K-l (11) 

and, 

*2*+l =d n~ a i X 2n + X ln+2 ) n = ~IA-^ ~ 2 (12) 

where a = -0.5 and (3 = 0.25. The multiplier and adder units of Fig. 1 will be modified 
accordingly to compute the pixels in accordance with equations (11) and (12). This filter 
takes as input the coefficients s 0 to s k -i and do to d k ., , and also requires as input the 
coefficients d./ and x. 2 from the previous band, which are stored in the partial register 109 
or external buffer 107. 

The aforementioned embodiment of Fig. 1 can also be modified without departing 
from the spirit of the invention to other wavelet filters such as a 9/7 wavelet filters. 

In this modified embodiment, the 9/7 forward transform wavelet filter computes the 
low pass coefficients and high pass coefficients in accordance with the following 
formulae: 

d' n =x Sn+1 + a(x s „+x 2n ^) n = 0,1,...,K-1 (14) 
So=* 2 * + P(d' n . 7 +d' n ) n = 0,1,...,K-1 (15) 
d n =d'„+ r (s„+s' n ^) n=-1,0,.~,K-2 (16) 
s tt =s' a +S(d n _ 1 +d n ) n = -1,0,...,K-2 (17) 
where d„ are all the high pass coefficients, s„ all the low pass coefficients, d' n and s '„ are 
intermediate values, a = -1.5861, p = -0.052980, y = 0.88291, and 5= 0.44351. The 
multiplier and adder units of Fig. 1 will be modified accordingly to compute the 
coefficients in accordance with equations (14) to (17). This filter takes as input pixels x, 
to xn and also requires as input the pixel x 0 , the intermediate values d'.j and s'.j, and the 
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high pass coefficient d.2 from the previous band, which are stored in the partial register 
109 or external buffer 107. 

In another modified embodiment, the 9/7 inverse transform wavelet filter computes 
the original pixels from the low pass coefficients and high pass coefficients in accordance 
with the following formulae: 

s' n =s n - 5(d nJl +d n ) n = 0,1,...,K-1 (18) 

d n =d n - y(s n +s' n+1 ) n~-1,0,...,K-2 (19) 

x 2 „=s' n - J3(d'„. y +d' n ) n = -1,0,...,K-2 (20) 
x sn *i=d'„-a(x 2lt +Xz fl .2) n = ~2,-1,...,K-3 (21) 
where d„ are all the high pass coefficients, s„ all the low pass coefficients, d'„ and s'„ are 
intermediate values, a = -1.5861, p* = -0.052980, y = 0.88291, and 5= 0.44351. The 
multiplier and adder units of Fig. 1 will be modified accordingly to compute the pixels in 
accordance with equations (1 8) to (21). This filter takes as input coefficients so to s^./ and 
do to dk-i and also requires as input the pixel x^, the intermediate values d'. 2 and s and 
the high pass coefficient d*j from the previous band which are stored in the partial 
register 109 or external buffer 107. 

The aforementioned embodiment shown in Fig, 3 can also be modified in similar 
fashion to other wavelet filters, such as a 5/3 inverse transform wavelet filter or 9/7 
wavelet filters. 

Industrial Applicability 

It is apparent from the above that the embodiment(s) of the invention are applicable 
to the computer and data processing industries. In particular, the embodiments have 
application in digital data compression of digital images and videos. 

The foregoing describes only one embodiment/some embodiments of the present 
invention, and modifications and/or changes can be made thereto without departing from 
the scope and spirit of the invention, the embodiment(s) being illustrative and not 
restrictive. 
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